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Abstract 



Differential cross sections for the electro-disintegration process e + 4 He — > 
3 H + p + e' are calculated, using a model in which the final state interaction 
is included by means of a nucleon-nucleus (3+1) potential constructed via 
Marchenko inversion. The required bound-state wave functions are calcu- 
lated within the integrodifferential equation approach (IDEA). In our model 
the important condition that the initial bound state and the final scattering 
state are orthogonal is fulfilled. The sensitivity of the cross section to the input 
p 3 H interaction in certain kinematical regions is investigated. The approach 
adopted could be useful in reactions involving few cluster systems where ef- 
fective interactions are not well known and exact methods are presently un- 
available. Although, our Plane-Wave Impulse Approximation results exhibit, 
similarly to other calculations, a dip in the five- fold differential cross-section 
around a missing momentum of ~450MeV/c, it is argued that this is an ar- 
tifact of the omission of re-scattering four-nucleon processes. 

PACS numbers: 21.45. +v, 25.10.-hs, 25.30.Fj 



I. INTRODUCTION 

In recent years, electro-disintegration processes in light nuclei have attracted much ex- 
perimental and theoretical attention. In the three-nucleon case, early calculations of the 
electron- proton coincidence cross section, i.e., of the scattered electron detected in coinci- 
dence with the ejected proton, were performed by Griffy and Oakes using analytical wave 
functions without inclusion of the final state interaction (FSI) . Since then, few-body tech- 
niques were used to obtain exact wave functions. Using the Faddeev formalism, Lehman 
calculated the two- and three-fragment electro-disintegration of 3 H and 3 He for a separable, 
Yamaguchi-type potential. In the final state the interactions between the ejected nucleon 
and the spectator pair were neglected. 
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Heimbach et al. || also studied the two-fragment electro-disintegration of 3 He and 3 H, 
but with the inclusion of the FSI. These calculations were based on a method introduced 
by Barbour and Phillips M, and extended by Gibson and Lehman ||, for the photo- 
disintegration of 3 He. The equations in this approach were derived by making use of the 
Alt-Grassberger-Sandhas (AGS) formalism for the three-nucleon system ||. As in Ref. [0, 
a spin-dependent, s-wave separable potential of Yamaguchi form was used in these calcu- 
lations. In contrast, van Meijgaard and Tjon |7| employed the semi-realistic Malfliet-Tjon 
(MT I-III) potential || represented in separable form via the unitary pole expansion (UPE). 

Realistic interactions have also been used to determine the electro-disintegration cross 
sections. Schiavilla || , for example, employed the Variational Monte Carlo ( VMC) method 
using the Urbana |L0| and Argonne [|ll|] nucleon-nucleon (NN) interactions for determining 



the bound state wave functions in the electro-disintegration of 3 He. His method involved 
constructing the final scattering states to be orthogonal to the initial bound state. In a 
completely different method, Laget |T2| employed a diagrammatic expansion of the scattering 
amplitude, which can be used to account for the most important multi-particle effects in 
the nuclear medium such as FSI and meson exchange currents (MEC). In a more rigorous 



approach Ishikawa et al. |13| calculated the electron-induced two-fragment disintegration 
of 3 He by using a 34-channel treatment based on Faddeev-type equations in momentum 



space. In their work they employed realistic two-nucleon interactions, namely the Paris [14 



and Bonn-B |T5| potentials. Later on, Golak at al. |16| extended the work of Ref. (T^J by 



including higher angular momentum components of the NN force (up to j = 2). 

In contrast to the three-nucleon electro-disintegration, calculations for the four-nucleon 
system are computationally much more demanding and to date no calculations have been 
performed to reliably include scattering states in a rigorous way. One therefore resorts to 
using the approximation 

<-><<!,,; tfA-xl- ^(qJi^l. (1) 

which is often employed in nuclear physics, and has previously been used in the photo- 
disintegration |L7f and electro-disintegration fl8| of 4 He. Here ^~- ) (q Q ; denotes the full 



scattering state associated with the channel state (q a \(^A-i\ in which the ejected proton 
is moving with momentum q a relative to the residual nucleus represented by the bound 
state (\Pyi_il. In such a model approach extreme care must be taken to ensure that the 
approximate scattering states in ([!]) are orthogonal to the initial bound state \^ a), 



(-) 



<q a |(*A-i|*A> = 0. (2) 



Non-fulfillment of this condition means that the effective scattering state, although asymp- 
totically correct, can be different from the "true" wave function in the interior. This in 
turn means that the off shell characteristics of the effective interaction employed can be 
completely wrong. This point is crucial since in photo- and electro-processes we are con- 
cerned with overlap integrals that include also the interior part of the wave-function. That 
condition is often simply ignored or some kind of an orthogonalization method is employed 
to ensure that it is fulfilled (see, for example, Ref. ||). The approximation ([J) together with 
the condition (|2|) has been used by Schiavilla [TEj within the VMC method (in a fashion 
similar to the one employed in 3 He electro-disintegration) using the Argonne two-nucleon 



and Urbana- VII three-nucleon |1| interaction models. The optical potential of van Oers et 
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al. pO| was employed in this approach to describe the 3+1 interaction. Laget J2T[] included 
the FSI and the meson exchange currents (MEC) using a diagrammatic expansion similar to 
the three-nucleon case. Nagorny et al. P^] included the electromagnetic field in the strongly 
interacting system in a relativistic and gauge invariant way, with the FSI being incorporated 
via the pole piece of the p 3 H^ p 3 H scattering t-matrix |23|. In a more recent work |2JJ 



the Schiavilla, the Laget, and the Nagorny methods were employed to study newly obtained 
data by van Leeuwe [53] . The results thus obtained were interpreted as an indication that 



the dip could be due to FSI effects and contributions from two-body currents. 

The main findings in the aforementioned investigations of three- and four-nucleon dis- 
integration processes were that (a) in certain kinematical regions there is an extreme sensi- 
tivity to the nuclear forces; (b) there are important differences between the results obtained 
when using different reaction models, and (c) in the four-nucleon case all model calculations 
within the Plane-Wave Impulse Approximation (PWIA) exhibit a dip in the five-fold dif- 
ferential cross-section around a missing momentum of ~450 MeV which, being absent from 
the experimental data, seems to be an artifact of the models. 

In this paper, we investigate the above sensitivities of the cross section to the input 
nucleon-nucleus interaction in the exclusive electro-disintegration process 

e + 4 He — ► 3 H + p + e' , (3) 

where the scattered electron and the ejected nucleon are assumed to be measured in coinci- 
dence. 

In our model the interaction between the outgoing proton and the residual nucleus is 
determined by using the nucleon-nucleus scattering data to generate an effective potential 



via the Marchenko inverse procedure p6|-p8|. The potential obtained would be unique, 



provided the phase shifts were known at all energies. Since theoretically the phase shifts can 
be calculated reliably only below the three-body break up threshold while experimentally 
they are available from a limited number of phase shift analyses, extrapolations to higher 
energies are required. The most obvious choice appears to be the one leading to states which 
fulfill the orthogonality condition given by Eq. (Q). Thus the dependence of the interaction 
on the extrapolated phase shifts and the resulting effect on the cross section is studied. Such 
an approach could be useful in reactions involving few cluster systems where the effective 
interaction needed to construct the scattering wave function, contains a lot of inherent 
ambiguities mainly stemming from the limited scattering data. In contrast, the bound 
states of the clusters involved can, nowadays, be calculated quite reliably using a variety 
of methods (Hyperspherical Harmonics, Variational etc.) with realistic or semirealistic NN 
forces. In our case the three-body and four-body bound-state wave functions are calculated 
in the integrodifferential equation approach (IDEA) P5|-|3"T|] using the MT I-III potential || . 

This paper is organised as follows: In Sect. [TT] we describe the cross section for the 
electro-disintegration 4 He in configuration space. In Sect. [I I J the electro-disintegration of 
4 He is discussed and the results of the model described in this paper are compared with 
the available experimental data and other theoretical results. Our concluding remarks are 
presented in Sect. |TV] 
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II. THE CROSS SECTION 



The electron-proton coincidence cross section is given by 
d 3 a a M Pf 



dEfdVt p dVt e 



(^c) 3 (2vr) 3 



n 

AEiEf cos 2 - 



\M(q c 



(4) 



where <7m is the Mott differential cross section, 



0"M 



4 2 

e cos 



AEf sin 4 - 
2 



(5) 



Ei(Ef) is the energy of the incoming (outgoing) electron and p/ is the relativistic density of 
states. The transition matrix is given by 



M(q c 



.(-) 



(6) 



where H is the effective Hamiltonian describing the interaction between an electron and the 
nucleons. The fragmentation considered here is of the 3+1-type while the ejected proton 
moves away with momentum q a with respect to the residual bound cluster described by the 
bound state \1/ a-i- 



The Hamiltonian for the interaction between an electron and A nucleons, is that of Mc 



Voy and van Hove |32j which has been previously employed in the electro-disintegration of 
the trinucleon system . This effective Hamiltonian, correct to order h 2 q 2 /M 2 c 2 , is 



H 



Aire' 



A* J 



.7 = 1 I 



[( Pj ■ a) e~ i( l- x * + e~^- x i ( Pj ■ a)] 



—i 



2M 

F 1N (q 2 ) + KF 2N {q : 



2\1 



2M 



<jj • (q X a) e 



iqxj 



8M 2 



[F lN {ql) + 2KF 2N {ql)} K) 



(7) 



where Xj and p- are vectors denoting the position and momentum of the j-th nucleon; <7j 
is the nucleon spin operator, a. is the electron's Dirac operator acting on the free electron 
spinors \m) and \uf), while q 2 is the exchanged four-momentum squared; Fin and F 2 n are 
the Dirac and Pauli form factors of the nucleon, k is the anomalous moment of the nucleon 
in nuclear magnetons, and M is the nucleon mass. For a proton knock-out the transition 
matrix, Eq. (J6|), is written 



M = (u f \(f\-(l + r) J 



^[F lp (ql)+^ p F 2p (ql)} 



e~ i( l- x ^ A )\ Ui ) 



A 



(u f \(f\cx 



+ -(1 + t) 3 z 



2 (1 + ^'2t 



2M 



cr 7 x q e 



-iqx 3 



\^>A)\Ui). 



(8) 



Here the subscript p refers to the proton. The transition matrix can be written in the 
convenient form as 



M = - 



u 



f\ui)Q - (uf\a\ui) ■ J 



(9) 



where 



Q = F! h {l + ql/SM 2 )(f 
J 
J el 



£e-^(l + T3 ) 



jel + jmag^ 
A 



f 



1 



£(^/2M)( Pj e-^ + e -«I-*ip.) (1 + T ,) 



(z/2M)F mag / 



E^x4(l 7- :{ ), 



(10) 



In these equations, T3 is the nucleon isospin operator, and F^ h and -F mag are the proton charge 
and magnetic form factors defined by 

F^ = F lp +{ql/AM 2 )K p F 2p 



mag 



(11) 

F lp + n p F 2p . (12) 

The analytical fit to the proton form factors F lp and F 2p , given by Janssens et al. is 
used in the calculations. 

Squaring the matrix element and summing and averaging over electron spins yields 

\ E \M\ 2 = [mE f + q l)QQ*-qlJ.r 

electron spins 

+2(kf ■ J)(ki ■ J*) + 2(kj ■ J*)(ki ■ J) 
-2E f {(k i -3)Q* + (k i -r)Q} 

-2E f {(k r J)Q* + (k r r)Q}. (13) 
Substituting this expansion into Eq. (f|) results in 



d 3 



a 



0"M 



\Pp\ E p 



dE f dtt p dQ e m C ) 3 (2vr) 3 , E p p p p A ^ 



l \Q\ 2 - X - sec 2 °- (Q*3 + 3*Q) ■ (h + k f ) 

+ - sec 2 - (J ■ fcjJ* ■ k f + J ■ k f 3* ■ hi) 
+ IJI 2 tan 2 - 1 . 



(14) 
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The evaluation of the coincidence cross section is thus reduced to the evaluation of the 
nuclear matrix elements Q and J, which depend on the choice of the initial and final wave 
functions. 

The inclusion of the FSI presupposes that the full solution of the scattering state be 
determined. Although this has been achieved in the three-nucleon case, it has not yet been 
accomplished in the four-nucleon case beyond the break up threshold. To circumvent this 
problem we employ the model outlined in the introduction in which the FSI is taken into 
account via an optical potential treatment of the relative motion of the outgoing clusters. 
In this particularly simplifying modification the plane wave (q a \ is replaced by scattering 
states ^ (q a \ generated by a p— 3 H potential obtained using the Marchenko inverse scattering 
procedure [p6|-|2~8[l- In other words, instead of solving a four-body integral equation providing 



us with the full scattering state ^(q a ; ips], we use the approximation given by Eq. (|T|). 

In implementing this model, we found it was numerically more stable to solve the integral 
equations in momentum space, in a similar fashion to the method employed by Fiedeldey 
et al. |35| in the photodisintegration of 4 He. In this method, the scattering state and the 
channel state can be projected onto each other by means of M0ller operators [Q. The 
resulting equation is of the Lippmann-Schwinger-type and can be written schematically as 

\M) = \B) + VG \M). (15) 

\B) is the plane-wave (Born) term, which can be evaluated in either configuration or mo- 
mentum space using the wave functions described below. 

The three- and four-nucleon bound state spatial wave functions are obtained using the 
IDEA. In this method, the A-nucleon wave function \l/(cc) is written as a sum of subampli- 
tudes *ff(x) = J2i<j<Ai J (. r iji r ) obeying the Faddeev-type equation 

[T + A{A ~ l) V Q {rmr ih r) = -[V (r^) - V (r)\ £ ^,r) (16) 

i<j<A 

where r is the hyperradius of the system, r 2 = -| T % and Vo(r) is the hypercentral potential 
which is the first term of the Potential Harmonic expansion of the interaction. The effects 
of the all important higher partial waves are approximately included via the use of the 
hypercentral potential (for more details of the IDEA see Refs. |29|-|3l|). 



III. RESULTS 

Reliable experimental phase shifts required in the construction of the effective potential 
for the 3+1 fragmentation of 4 He are not available. Furthermore, the energy region in which 
they are available has little or no overlap with the region needed for our calculations. On the 
theoretical front, phase shifts beyond the three-body break up have not yet been obtained. 
Therefore, in implementing our model we employed phase shifts obtained from the optical 
potential used by Schiavilla fli8| , |2"0|| . Since the phase shifts produced by this potential become 



negative for energies greater than 200 MeV, we had to extrapolate them to higher energies. 
In Fig. [I| four such extrapolations of the real phase shifts to higher energies are plotted. The 
corresponding real parts of the Marchenko potentials that describe the p+ 3 H interaction are 
plotted in Fig. 0. The local potentials produced by the Marchenko procedure are complex 
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and unique for each set of data. They are phase equivalent up to approximately 300 MeV, 
producing phase shifts with slightly different high energy behavior. The extrapolation Ext. 
2, produced an overlap integral that was, for all practical purposes, zero. 

The various kinematics used are given in Table |. In Fig. [3] we have plotted our cross 
section results for the kinematics R together with the experimental data of Ref . [f25| . These 
results are very close to those obtained by Laget ||21|| . It is obvious that the cross section 
is only weakly dependent on the nucleon-nucleus interaction, all extrapolations yielding 
virtually identical results. 

The effect of the high-energy behavior of the phase shifts on the cross section and the 
importance of the orthogonality condition given by Eq. (Q), is more clearly demonstrated 
in the plots of the cross sections for the kinematics (A,B,C,D) and kinematics S, shown in 
Figs. ^ to |6]. The PWIA prediction for kinematics (A,B,C,D) shows a minimum, in fact a 
zero, in agreement with the calculations of Laget and Schiavilla. At missing momenta below 
300MeV/c, all calculations show good agreement with the data, but in the region of the 
dip the calculations underestimate the data considerably. Beyond missing momenta of ~ 
550MeV/c the agreement of the calculations with the data is again fair. The best result 
is the one obtained for Ext. 2 which fulfills the orthogonality condition (Q). In Fig. |5|, we 
compare the cross section produced using this extrapolation with the results of Laget [ 2~I | ; 
it is seen that the agreement is good. 

The PWIA cross section calculated for kinematics S, shown in Fig. |6], also exhibits a dip 
at ~450MeV/c. Including the FSI causes a partial filling of this dip. The position of the 
minimum and the value of the cross section in this region once again depends strongly on 
the input nucleon-nucleus and the NN potentials. 



IV. CONCLUSIONS 

The results described above clearly indicate a sensitivity of the cross section in certain 
kinematical regions to the input nucleon-nucleus interaction and to the input NN potential. 
In our model, the various nucleon-nucleus (3+1) interactions employed for this investigation 
were phase equivalent up to approximately 300 MeV. They all reproduced the correct binding 
energy for the four-nucleon bound state. However, beyond 300 MeV, the potentials produced 
phase shifts that differed slightly. The corresponding nucleon-nucleus potentials provided 
scattering states that were not always orthogonal to the initial bound state. This in turn 
is translated into large differences in the cross sections in certain kinematical regions. This 
result was not surprising, as already in 1970, Fiedeldey |38j investigated the dependence 
of the triton binding energy on the high-energy part of the phase shift. He concluded 
that arbitrary variations of the NN phase shift at high energies (E\ a ^ > 300 MeV) can 
produce large differences in the triton binding energy. This sensitivity is also manifested 
in the electro-disintegration cross sections which implies that we need a much clearer idea 
of what constitutes a physically acceptable extension of the phase shift to higher energies. 
This sounds a warning that in model calculations care must be taken when using effective 
interactions, the most important constraint being the condition that the initial bound state 
and the final scattering state be orthogonal. 

The dip around 450MeV/c, which is present in all model calculations, is something of a 
surprise. At missing momenta less than 300MeV/c all calculations show a good agreement 
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with the experimental data, i.e., the PWIA performs reasonably well in a region where the 
FSI could be expected to be more important. At missing momenta beyond the dip area, 
where the MEC become important, the agreement with the data is fair. Thus the zero in 
the PWIA cross-section is not necessarily a manifestation of strong FSI effects or MEC. 
It can simply be attributed to the vanishing Born term which contains only genuine 3 + 1 
components. The Born term should, however, be treated in a more rigorous way using the 
AGS formalism which allows for the effects coming from the 2+2 channel. The inclusion of 
this channel should result in an improvement of the the bound state wave functions and of 
the Born term. This in turn is expected to remove the apparent zero in the PWIA cross 
section. Investigations in this regard are currently under way. 
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TABLES 



TABLE I. Kinematics R and (A,B,C,D) from Ref. [25], and kinematics S from Ref. [37]. 





Ei [MeV] 


E f [MeV] 


e 


q [MeV/c] 


R 


524.9 


423.9 


49.60° 


408 


(A,B,C,D) 


525.0 


310.0 


49.60° 


401 


S 


560.0 


360.0 


25.00° 


278 
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FIG. 3. The two-fragment electrodisintegration cross section of 4 He as a function of the 
missing momentum for kinematics R of Table 1. 

1 



a 

■d 

a. 

a 



"i 1 1 1 r~ — r — rr~ i 

Experiment [25] i — i — i 

PWIA [21] 

Ext 2 . 

5"<hL Ext 3 

Ext 4 




250 300 350 400 450 500 550 600 650 
p m [MeV/c] 



FIG. 4. The two-fragment electrodisintegration cross section of 4 He as a function of the 
missing momentum for kinematics (A,B ; C,D) of Table 1. 
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FIG. 5. A comparison of our two-fragment electrodisintegration cross section of 4 He with that 
of Laget as a function of the missing momentum for kinematics (A,B,C,D). 
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FIG. 6. The two-fragment electrodisintegration cross section of 4 He as a function of the 
missing momentum for kinematics S of Table 1. 
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